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EFFECT ON METEOR FLIGHT OF COOLING BY RADIATION 


AND ABLATION 

By H. Julian Allen, Barrett S. Baldwin, Jr., 
and Nataline A. James 
Ames Research Center 


SUMMARY 

gjc A 

The parameters of meteor flight are evaluated taking account of surface 
radiation as well as vapor ablation, Results are presented for four meteoroid 
materials: tektite, stone, pumice, and iron. The calculations are based on a 
spherical shape with a range in radii of 10” to 5x10 ” 3 meters. Entry veloci- 
ties from 11 to 70 kilometers per second are considered. Because of an approx- 
imate similarity a large range of entry angles is covered. The validity of 
this approximation and of a number of assumptions is investigated. A signifi- 
cant result of the investigation is the rediscovery of surface radiation as a 
possible explanation, instead of fragmentation, for 
the light curves of faint meteors. 

INTRODUCTION 



It is usual in meteor physics to assume (e.g., see refs. 1 and 2) that 
the fieat energy resulting from the collision of air molecules with the surface 
of a meteoroid is completely accounted for by vapor ablation of the surface. 
This assumption is a reasonable one for meteoroids of centimeter size and 
larger. Whipple (ref. 3) showed that if the size of the body at atmosphere 
entry is sufficiently small, say of the order of a micron or less, then the 
reradiation of heat from the meteoroid alone will suffice to keep the surface 
cool enough that essentially no ablation will occur. It is the purpose of 
this paper to analyze the motion and heating of meteoroids of sizes between 
the two extremes in order to answer, for a wide range of speeds, such questions 
as: (l) At what size does melting occur and over what range of altitudes does 

the molten state exist? (2) Under what conditions will a nonrotating body 
with a molten front face have a solid rear face? (3) What is the maximum size 
for which one may expect free molecular flow to exist? (4) Up to what size 
will a molten drop remain essentially spherical? ( 5 ) When will the dynamic 
pressure overbalance the surface tension and break up a molten drop? (6) For 
what meteoroid materials will vapor ejected from the front face significantly 
alter the drag force? (7) How does the size of the meteorite vary with the 
initial size of the meteoroid, and, in particular, for irons, is the calculated 
meteorite size -abundance relation in agreement with that obtained from cosmic 
dust collections from polar ice caps? (8) At what altitude does the speed fall 
to a trivially low value? and finally, (9) What is the behavior of the luminos- 
ity of these bodies with respect to both maximum magnitude, and magnitude var- 
iation with altitude? 



In the analysis four meteoroid materials have been considered. In 
addition to the usual stone and iron, we have included tektite, which Chapman 
(ref. 4) has indicated to be of lunar origin, and a low density material with 
the character of pumice. Much of the needed information on radiative and 
vaporization properties of meteoroid materials is contained in reference 5* 

The important physical characteristics assumed for these materials are listed 
in table I. The symbols in the analysis are those used in meteor theory and 
familiar to astrophysicists; the system of units is the meter -kilogram-second. 


SYMBOLS 


The quantities used in the analysis along with the units in which they 
are expressed are the following: 

a meteoroid frontal area at any time t, corresponding to altitude h, m 2 

a dimensionless meteoroid frontal area, a/a 

a^ meteoroid frontal area at atmosphere entry, m 2 

b sticking coefficient 

c 1} c 2 constants in the vapor pres sure -temperature equation 
d diameter of vapor molecule, assumed to be 2xl0 -1 ° m 

h instantaneous altitude, m 

I luminous intensity, t q kg m 3 sec -4 

I max maximum luminous intensity, t q kg m 3 sec -4 
i luminous intensity coefficient, sec -1 

i maY maximum luminous intensity coefficient, sec -1 

Kn Knuds en number, ij r 

k thermal conductivity of the meteoroid, W m -1 °K -1 

L Avogadro's number, 6.02252xl0 23 mole -1 

l mean free path of air at altitude h, m 

M molecular weight of meteor vapor, kg mole -1 

m meteoroid mass at time t, corresponding to altitude h, kg 

m dimensionless meteoroid mass, m/m 
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meteoroid mass at atmosphere entry, kg 
vapor mass from front face of meteoroid, kg 

number of vapor molecules from the front surface of the meteoroid 
flowing through the stream tube of area a 

meteoroid saturation vapor pressure at temperature T, N m" 2 

front surface meteoroid vapor pressure at temperature T x , N m" 2 

rear surface meteoroid vapor pressure at temperature T 2 , N m" 2 

heat input in kinetic energy units, J 

gas constant (8.31^3/M), m 2 sec" 2 °K -1 

meteoroid radius at any time t, corresponding to altitude h, m 

dimensionless radius, r/r 

final (meteorite) radius at sea level, m 

meteoroid radius at atmosphere entry, m 

surface area of the meteoroid at any time t, m 2 

temperature, °K 

melting temperature of meteoroid, °K 
front surface temperature, °K 
rear surface temperature, °K 
time, sec 

meteoroid velocity at any time t, corresponding to altitude h, 
m sec" 1 

dimensionless meteoroid velocity, V/V 
meteoroid velocity at atmosphere entry, m sec" 1 
Weber number, pV 2 r/r s 

distance measured ahead of the body in the direction of motion, m 
zenith angle (trajectory angle measured from vertical), rad 
inverse of over-all scale height for density, m" 1 



P’ inverse of local scale height for density, m -1 

7 drag coefficient 

7 v coefficient of the drag increment due to meteoroid vaporization 

e radiation emissivity (assumed to be 0.92) 

£ energy required to vaporize a unit mass of the meteoroid, m 2 sec -2 

(J kg' 1 ) 

q kinetic to ablation energy ratio, AVe 2 /l27£ 

A heat -transfer coefficient 

| dimensionless air density factor 

p air density at altitude h, kg m -3 

p dimensionless air density, p/p Q 

P m meteoroid density, kg m“ 3 

p Q sea-level air density, 1.225 kg m -3 

p v ^ front surface meteor saturation vapor density, kg m" 3 

a Stefan -Boltzmann constant, 5-6697x10 8 W m -2 °K 4 

r 0 luminous efficiency factor, I sec 4 kg' 1 m" 3 

T g liquid meteoroid surface tension, N m -1 

\|r molecular screening factor (ratio of total area of vapor molecules in 

the stream tube ahead of the meteoroid to the meteoroid frontal area) 

cj mass evaporation rate per unit area, kg m -2 sec -1 

cj-l front surface evaporation rate per unit area, kg m -2 sec' 1 

co 2 rear surface evaporation rate per unit area, kg m -2 sec' 1 

ASSUMPTIONS 

(1) The weight can be ignored and the earth assumed flat in the motion 
calculation. 

(2) For the motion and ablation calculations, the meteoroid is initially 

a nonrotating sphere before entering the atmosphere and it remains a nonrotat- 

ing sphere as it ablates. 
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(3) The flow is free molecular. The limits for which this assumption is 
applicable are examined. 

(4) Change in the heat stored within the body can be ignored compared 
with the heat removed by radiation and by ablation of the surface as vapor. 

(5) A temperature difference between the front and rear surfaces of the 
body may exist. For purposes of calculation of this effect, the body is 
treated as a slab of thickness equal to the radius. 

(6) The variation of atmospheric density and mean free path with altitude 
is given by the 1959 U. S. Standard Atmosphere (ref. 6). 

Most of these assumptions are usual in meteor theory. The justifications 
for them have been so fully discussed in the literature that no additional 
remarks are felt to be needed here. 


ANALYSIS 


Since, for the most part, the analytical development herein is that of 
the classical theory of meteors (see, e.g., ref. 5 )* derivations of formulas, 
in the main, are included only where they are new or obscure in the literature. 

In accordance with the first assumption the meteor flight path is a 
straight line, so that time can be expressed in terms of altitude through 


dt 


dh 

V cos Zr 


The equation of motion for the meteoroid is, with assumption (l), 


( 1 ) 


With assumption (2) 


and equation 


m ^ = -7pV 2 a 
(l) this becomes 


av = 3 p 0 7 / pv\ 

^ ^P m r 00 cos 2r 


The heat input rate is given by 


dq 

dt 


= J ApV 3 a 
2 


In accordance with assumption (4) the heat balance gives 

= k ecrs(T 1 4 + T 2 4 ) + J SsK + u 2 ) 


( 2 ) 

(3) 

W 

(5) 
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(6) 


With assumption (2) this becomes 

(5V 3 ) . €a(T x 4 + T 2 4 ) + CK + “ s ) 

The heat received by the rear face must be conducted from the front face so 
that in accordance with assumption ( 5 ) 



€CT 2 4 


+ f> 2 


( 7 ) 


The mass loss by ablation is 

~ = -V cos Z R g = - \ sK + ui 2 ) (8) 

With assumption (2) this becomes 


311 " 2 VU» 308 Z R 


( 9 ) 


To use equations (6), ( 7 ), and (8) it is necessary to express the evaporation 
rate in terms of the surface temperature. For a Maxwellian distribution of 
efflux velocities (see ref. p. 2b) 


and 


bPvi 

«i = r - ... 

•7 2jtRT x 


Cl+( C 2/ T l) 
bxlO 

<s / 2?t RT^ 


( 1)2 


^Pv 2 
•J 2rtRT 2 


bxlO 


Ci+(c 2 /T 2 ) 


n/ 2jtRT 2 


(10a) 


(10b) 


The quantity b is the so-called sticking coefficient (see ref. 7> P« 85 and 
ref. 8). For iron b is near 1.0. For solid stone b can be small, but 
according to reference 8, its value should be near 1.0 for liquid stone. 

Since evaporation of stone from the solid state is not appreciable we shall 
set b = 1.0 throughout the present investigation. 

To use equations ( 3 ) and ( 5 ) we must also provide the heat transfer and 
drag coefficients. With assumption ( 3 ) and complete accommodation of molecu- 
lar impact of the air on the surface 

A = 1 (11) 

The drag coefficient would also be unity for complete accommodation but 
when the vaporization of the front surface differs from that of the rear sur- 
face, then a drag increment due to vaporization will be present in the amount 

7vPV 2 a = | a(p Vi - p Vg ) (12) 
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or , 


7 y = 


1 


2PoV 


" 10 Ci + (c2/Ti) _ 10 Ci+(c 2 /T 2 )- 


(13) 


and the total drag coefficient will thus he 

7 = 1 + 7 V (l 1 ^) 

Using the foregoing equations one can calculate the entry characteristics 
for a meteoroid in free molecular flow. For a nonahlating body, the free 
molecular flow limit occurs when the body size at any altitude becomes compar- 
able with the mean free path of the air molecules at that altitude. In other 
words the Knudsen number 

Kn = i (15) 

must be substantially larger than unity (say, about 5). 


The free molecular flow limit for an ablating body may also be due to 
ablating vapor molecules colliding with a sufficient fraction of the air mole- 
cules ahead of the body, thus, effectively shielding the body from these air 
collisions. The screening factor, which is a measure of the departure from 
free molecular flow, may be found as follows: The vapor density at the front 

surface of the body is 


10 G1+(C2 / Tl) 

Pv i " RTi “ RT X 


(16) 


and at any time the total mass of vapor in the stream tube of unit area ahead 
of the body is 


m Yl 


a 



The number of vapor molecules per unit area is then 




c x +(c 2 /T x ) 
LrlO _____ 

RT X 


(17) 


(18) 


The diameter of each molecule is d, so that the total cross section of 
all vapor molecules in the stream tube of ■unit area is 




= rtd 



rtd Lr, 


CO. 


R 


Cx+(c 2 /T 1 )- 1 


rlO 


(19) 
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If the free molecular flow assumptions are to apply, the screening factor, 
\|/, should be small compared to unity (say, 0.1 or less). However, it may well 
be that even for much larger values, the departure of the heat transfer and 
the drag coefficient from free molecular flow values may not be great. In the 
results presented later there are plots showing the values of entry parameters 
that lead to maximum values of \|r cos Zr equal to 0.1, 0.3^ 0.5, and 1.0. A 
similar criterion for screening due to vapor molecules is given in reference 7 
on pages 35-39 and 82-85. 


A second limitation (see ref. 5) which applies to the analysis is that 
the meteoroid should stay essentially spherical during entry and must remain 
intact. When the dynamic pressure, p V 2 , becomes comparable with the internal 
pressure in the molten drop due to surface tension, the body becomes distorted 
from a spherical shape. Further increase of dynamic pressure flattens the 
drop, and when this flattening becomes too great, the drop becomes unstable to 
small disturbances and breaks apart. The Weber number relates the dynamic 
pressure and surface tension. 


We 


pV 2 r 



( 20 ) 


Chapman and Larson (ref. 9) note that for Weber numbers up to about 1.1, 
the departure from spherical shape should be trivial. The appropriate value 
of the Weber number for breakup is uncertain. Chapman and Larson used 

We = 2.6 (21) 

as the value for incipient disruption of a tektite drop, but it certainly is 
possible that a larger value would be more appropriate. The viscosity of the 
molten meteoric material plays an important role in this regard. For example, 
for molten iron the appropriate value would be perhaps 2 or less since the 
kinematic viscosity is very low. In the results presented later there are 
plots showing the values of the entry parameters that lead, to values of 
We cos Zr equal to 1.1, 2.5, and 4.0. 

Finally, the luminous intensity is defined by the usual relation (e.g., 
ref. 2) 



( 22 ) 


Let us define a luminous intensity coefficient as 


1 = 


21 


dm 


r o V °o 3 %> dt 

With assumptions (l) and (2) this becomes 


(23) 


i = 3V ro cos Zr ( vV 


(24) 


8 


RESULTS 


i 


The foregoing equations were solved on an electronic computer. Results 
I were obtained for the four meteoroid materials with the product r^ cos Zr 
varying from 0.01 p to 5 mm. 

The results are generally presented as if r^ cos Zr were, in fact, a 
proper similarity variable, which it is not when the temperature difference 
between the front and rear face becomes important, since this difference is 
determined by r rather than r cos Zb. A few check calculations indicated 
| that the errors introduced by this treatment were small, as will be shown 
later, so that for practical purposes the variable may be considered a proper 
one even though the results rigorously apply only when cos Zr is equal to 
unity . 

For all of the materials except tektite the assumed entry velocities range 
from 11 to 70 km/sec. For the tektites the range chosen was 11 to 30 km/sec. 
The higher speed results for stone and iron are only included to make the 
trends with speed self-evident. 

The results presented in figures 1 through 8 are for tektite, figures 9 
through 23 for stone, figures 24 through 32 for pumice, and figures 33 through 
41 for iron. For each of these materials at the various speeds and for the 
full range of r M cos Zr the following results are given: (l) front surface 

maximum temperature, (2) rear surface maximum temperature, (3) altitude range 
for which the molten state exists, (4) ratio of final size to size at entry, 
(5) altitude at which the speed has dropped to 500 m/sec, (6) values of 
j-max/ cos Zr, (7) altitude for maximum luminous intensity, and (8) the solid- 
liquid boundary as well as information on cos Zr and We cos Zr to indicate 
the likely limits of applicability of the results. Some additional figures 
are included which give other quantities that are useful to the discussion. 

In the case of luminous intensity, it is of value to compare the results 
of the computations with the predictions of a simplified theory which assumes 
that no radiation from the body occurs. The "ablation -only" theory in ques- 
tion is treated in the appendix. 


DISCUSSION AND CONCLUDING REMARKS 


In the following discussion of results we will, in a number of instances, 
compare the present machine calculated results with the predictions of usual 
meteor theory, which assumes that ablation alone is the cooling mechanism and 
that the atmospheric density varies exponentially with altitude: 


P = 



e- P h 


(25) 
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The inverse scale height, p, is not constant, in fact, hut varies as . 
shown in figure 42, which was plotted using data from the 1962 U. S. Standard 
Atmosphere (ref. 6). Because p does vary with altitude, it is to he 
expected that better agreement between the machine calculated results and the- 
ory would he obtained if in the latter at any given altitude of interest a 
local value of 3 (herein denoted as p*) appropriate to that altitude were 
used in the comparison. The value (3’ can he found as follows: Differenti- 

ating equation (25) with 3 as a function of h gives 


so that 


dp 

dh 


= -e 


•3h 


. dp 
+ h ~ 
dh 


= -P 3 + 


dp 

h dh 


3’ = 3 + 


, <13 


dp/dh 


(26) 


(27) 


Values of p* have been determined for the U. S. Standard Atmosphere 
1962 using the tabulated values (ref. 6) of p and dp/dh averaged over a 1 km 
height. These values are also plotted in figure 42. Note that the p* vari- 
ation with altitude is somewhat erratic even though the method used for obtain- 
ing dp/dh tends to smooth out the variation. That is, the densities given 
in this standard atmosphere do not fall on a smooth curve, and as we shall see, 
this fact reflects itself in certain of the results presented. 


When comparisons are made with theory assuming an average value of P, 
this value has been chosen as 1. 45x10 “ 4 m -1 . 

At the outset it is necessary to establish how well the product r^ cos Zp 
fits as a fundamental similarity variable. As noted earlier when front and 
rear temperatures differ, the product r M cos Zp is not a proper variable. 

The adequacy of treating it as such can be determined by examining the effect 
of changing cos Zr from unity to some small value for a material for which 
large differences between front and rear surface temperatures occur. Of the 
four materials studied, pumice is such a material. To find the results for 
pumice values of cos Zp of 0. 3 as well as unity were used. The deviation of 
results was significant only for those quantities shown in figures 25, 27, and 
32. The dashed curves show the effect of reducing cos Zp to 0.3 at 11, 30> 
and 70 km/sec. The differences are significant only for the ratio of maximum- 
to -minimum drag coefficient at the higher speeds and for the rear surface 
temperatures at the lower speeds. 

The numerical results indicated, that the effect of ablation on drag coef- 
ficient was important only for pumice and then only for the larger sizes. For 
that reason the ratio of maximum -to -minimum drag coefficient as a function of 
entry size and speed for the pumice only is shown (fig. 32). 


Consider, now, the variation of the maximum front and rear surface tem- 
peratures with entry size and speed (see figs. 1, 2, 9, 10, 24, 25, 33> and 
34). These variations are about what one would expect. The plot of figure 43 
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serves to assist in comparing the differences that do exist for the four 
materials at a common entry speed (30 km/sec). For sizes of 1 micron and less, 
radiative cooling predominates so that 

M i ApV 3 A = 4ecrT 4 A (28) 

dt 2 

Then it may he shown (e.g., using the analysis of ref. 10 which assumes no 
ablation) that for a given entry speed 

T Jax ~ PPm r oo cos % (29) 

This agrees with the results presented in figure 4-3 if the local inverse scale 
height, 3', is used in place of (3* 

As the sizes are increased into the range wherein ablation affects the 
results, it will be noted that the tektite temperatures become relatively high 
when compared with those of the other materials. This occurs because the 
ablation rate at a given temperature is less than that for the other materials 
(see values of c x and c 2 of table I). 

A comparison for each of the materials of the size at which the maximum 
front and rear temperatures depart from each other is somewhat surprising in 
that, although the thermal conductivities of the materials are greatly differ- 
ent, the "departure" size differences are not great. The conductivity of the 
iron is 400 times greater for the pumice but the departure size is only one 
order of magnitude larger. 

Before leaving the subject of maximum front surface temperature it should 
be pointed out that, except for stone and tektite, the results may be subject 
to future readjustment. For pumice the ablation parameters c x and c 2 may 
differ considerably from the values used (the same as for stone). For iron no 
allowance was made for heat addition resulting from oxidation, that theory 
indicates (ref. ll) and experiment shows (ref. 12) will occur. Thus for the 
larger sized irons the actual temperatures should exceed somewhat the values 
shown. 

One final comment concerning maximum rear surface temperatures is in 
order. For pumice only there exists a rather large range of entry sizes and 
speeds for which no rear surface melting occurs, although for tektite at the 
lowest speeds (which Chapman's analyses indicate to be the most likely) the 
rear surface would be a very viscous fluid for a narrow range of size. 

Examination of the plots of the ratio of final radius to initial radius 
as a function of entry size and speed (figs. 4, 12, 27, and 36) shows that the 
transition from the case for which radiative cooling dominates (vo/r^ = i) to 
that for which ablative cooling dominates takes place within a rather small 
change in entry size. This is particularly evident at the highest entry 
speeds wherein the transition is essentially_complete with one order of magni- 
tude change in entry radius . The plots of V versus r for stone at 30 km / sec 
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(fig. 17) are characteristic of the manner in which the size of all the 
materials changes. The dashed curve gives the theoretical variation when abla- 
tive cooling only occurs. (See appendix, eq. (A3).) 

The variation of the maximum luminous intensity coefficient with size and 
speed has some interesting features. (See figs. 6, 14, 29, and 380 As entry 
size decreases from the maximum of the study, the coefficient generally 
increases somewhat erratically to a maximum and then rather abruptly falls to 
a trivial value. This erratic behavior is evidently due to the fact that the 
standard atmospheric density variations with altitude given by reference 6 are 
themselves erratic, as was mentioned earlier. The peaking characteristic, 
which was not expected, cannot be explained by variations in the local values 
of p as is shown (fig. 18) for the example case of stones entering the atmos- 
phere at 30 km/sec. (On this and some subsequent figures the calculated 
points are included to show how the curves were faired.) The peaking behavior 
can be explained as follows: From equation (23) 

* - * I (30) 

and for ablative cooling only ( see appendix, eq. (Al) ) 

# - -tsvV ( 31) 

dt 

so that 

i - JffV (32) 

Bodies with entry radii of order 10" 5 m, which at the higher altitudes are 
primarily cooled hy radiation, remain larger throughout the trajectory than 
they would if ablation were the sole cooling means, as noted earlier. In par- 
ticular, the size at the altitude for maximum luminous intensity is larger 
(see fig. 19) hut, more importantly, the altitude itself (see fig. 20) is so 
much lower than predicted for the case of ablation cooling only that the cor- 
responding air density is considerably greater. Thus, although the velocity 
at peak luminous intensity is somewhat lower (see fig. 21), the product pV 6 !* 2 
is much larger than theory for ablative cooling only predicts. When the theo- 
retical values of luminous intensity using the local (3 are corrected for the 
change in this product, the results (fig. 18) provide some explanation for the 
peaking trend observed. Altogether, the peaking can be attributed to the 
existence of radiative cooling, which causes a deeper penetration into the 
atmosphere than would otherwise occur. 

The variations of maximum luminous intensity discussed above are only of 
academic interest since meteor magnitude is determined by the logarithm of 
intensity and the variation of the logarithms of the maximum luminous inten- 
sity coefficient is small as may be seen, for example, in figure 22. This 
figure shows how photographic or visual meteor absolute magnitude varies with 
altitude for the various sizes of stone bodies for the common entry speed of 
30 km/sec. A more instructive plot is that of figure 23 wherein the 
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magnitudes and altitudes are normalized with respect to the maximum mag nitude 
and altitude at maximum magnitude. On this plot one can compare the results 
with theory for ablative cooling only. The differences shown are about what 
one would expect. 

The appearance of figure 23 is of interest in another connection. Devia- 
tions of the light curve from the results of pure ablation theory of the type 
shown are commonly observed in the case of faint meteors (see refs. 13-15)* 

Cook (ref. 13) attributed such deviations to surface radiation. Later Jacchia 
(ref. 14) showed that fragmentation would cause a similar behavior. The frag- 
mentation explanation seems to have been widely accepted (see, e.g., refs. 5 
and 15). It is probable that the explanation based on surface radiation is 
the correct one for many of the faint meteors that have been observed. 

Some comments are in order on the limits of applicability of the analysis 
with respect to both departure from free molecular flow and freedom from drop 
breakup resulting from excessive dynamic pressure. 

Departure from free molecular flow occurs only as a result of interaction 
between vapor molecules from the front face with the air molecules since in no 
case did the Knudsen number fall below unity until that portion of the trajec- 
tory was reached wherein the speeds had dropped to very low values. The 
screening factor values, 4^^ cos Zr, are given in figures 8, 16, 31, and 40 
for the four materials, and they show that the analysis should be valid for 
cos Zr to about 100 p. 

Drop breakup does not impose a limit until the size is increased to 
nearly 1 mm as shown by the plotted Weber numbers. In fact, for slow speed 
tektites, and particularly for the pumices, drop breakup may not be a factor 
if, as indicated, the rear face of the body is solid or essentially so. 

A final comment is in order concerning a question raised in the introduc- 
tion regarding the compatibility of the calculated final size of irons pre- 
sented with the known iron meteorite size -abundance relation found from cosmic 
dust collection. 

Figure 4l presents the meteorite size (r Q cos Zr) as a function of entry 
size for entry speeds of 20 and 30 km/sec which should bracket the probable 
entry speed range for such meteoroids . As was noted previously the values of 
r 0 cos Z^ for the larger entry sizes might well be too big because in the 
analysis no account was taken of additional ablation resulting from heat addi- 
tion due to oxidation. With this restriction in mind and the fact that the 
average value for the cosine Zr can be shown to be two -thirds, the expecta- 
tion from these results is that the maximum radius of iron cosmic dust should 
be greater than a few microns and not much more than the order of 20 p. This 
result is in good agreement with the results of collection counts (see ref. 16). 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., March 18, 1965 
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APPENDIX A 


ANALYSIS OF ENTRY INTO AN ISOTHERMAL ATMOSPHERE 
WITH ABLATIVE COOLING ONLY 


If the meteoric body is cooled only by ablation, then 


C = - % *pV=>a 

dt 2 


(Al) 


Assuming an isothermal atmosphere together with assumptions (l) and (2) leads 
to 


dn 


PnH 


o' ’ 00 


(V 2 e' ph ) 


dh ep^r^ cos Z R 
With equation ( 3 ) then one obtains Spik* s equation 

_ -n(i-v 2 ) 

r = e 

where 


(A2) 


AV 


H = 


12r^ 


(A3) 


(A4) 


From equations ( 3 ) and (A3) one can show that 


U = e 


-» h . °° s Z R { 8 -n[g (l ) . iTW 2 )]} 


37P r 


(A5) 


where Ei denotes the exponential integral (for tabulated values see ref. 17 ). 
Substituting equations (A2), (A3), and (A5) in equation (24) gives 


- 3 tj 2 qV 2 , -2.3 

3 e e ) Ft?-? ( ^ ^ _ ir.-i ( 


Ei(rj) - Ei(qV 2 )] (A6) 


COS Z R tj* 

The maximum value of i/pV^ cos and the corresponding values of V, r, and 


I = 


7 P 


3p m r cc cos Z R 3 
have been calculated and are given in table II. 


= | e _ri [Ei(q) - Ei(qV 2 )] 


(A7) 
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Quantity 


Iron 


Stone 


Tektite 


Pumice 


M 

56 

0 

LT\ 

42 

50 

+ Cl 

9.607 

9.600 

13.301 

9.600 

+ c 2 

-1.6l2xl0 4 

-1.350xl0 4 

-2.5102xl0 4 

-1.350x10' 

m 2 sec" 2 

8.0ib<10 s 

7.98xl0 6 

14 . 98 x 10 s 

7.98x10 s 

P m > kg m" 3 

7 . 8 xl 0 3 

3 . 4 xl 0 3 

2 . 4 xl 0 3 

4 xl 0 2 

■ s , kg sec" 2 

1.20 

O.36 

0.36 

O.36 

W m' 1 °K” 1 

40.0 

2.0 

1.6 

0.1 

T f , °K 

1.8xl0 3 

I.8xl0 3 

1.8xl0 3 

l.8xl0 3 














AND RELATED QUANTITIES 


TABLE II.- VALUES OF gj, 

[Note 0.4 - 2 = 0.4xl0~ 2 ] 




X* * 

^-max 

V, 

i-max 

Imax/^Xx) cos Zp 

£ . 

^•max 
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1.0000 

0.84645 
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0.2222-0 

.2 

• 9500 

.86249 

.06422 
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• 5 

.8950 

.88217 

.13411 

.1579-0 

1.0 

.8361 

.90608 

.20897 
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. 7999 

.92260 

.25568 

.9598-1 

2.0 

.7760 

.93444 

.28728 

.7962-1 
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.7590 

.94323 

.30997 

.6792-1 

3.0 

• 7465 

.95002 

.32701 

.5912-1 

4.0 

.7294 

•95975 

.35084 

.4687-1 

5.0 

.7182 

.96634 

• 36668 

.3879-1 
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• 7104 

.97109 

•37795 

.3307-1 
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• 7046 

•97767 

.38638 

.2881-1 
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•97778 
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.98156 

. 40239 

.2076-1 
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Figure 1.- Front face maximum temperature, T x , for tektite 









Figure k.- Ratio of final to entry radius, ro/r^, for tektite. 





Figure 5 .- Final altitude (V = 500 m/sec) for tektite. 





Figure 7.- Altitude at maximum luminous intensity for tektite. 








Figure 9.- Front face maximum temperature, for stone 







= 70 km/sec 



Figure 11.- Front face melting zones for stone. 



Figure 12.- Ratio of final to entry radius, Tq/t^, for stone. 





= 70 km/sec 



Figure 13 . - Final altitude (V = 500 m/sec) for stone. 











= 70 km /sec 



Figure 15.- Altitude at maximum luminous intensity for stone. 



Figure l6.- Liin.i'ts of* validity of* analysis lor stone. 








Figure 18.- Maximum luminous intensity coefficient for stone at 30 km/ 





Figure 19-- Values of r at maximum luminous intensity for stone at 30 km/sec. 





30 km /sec /Theory for ablation cooling only 



Figure 21.- Values of V at maximum luminous intensity for stone at 30 km/sec 








Results for cos Z 





Voo = 70 km/sec 



Figure 26.- Front face melting zone for pumice 




Figure 2J.- Ratio of final to entry radius, Tq/y^, for pumice 





Figure 28.- Final altitude (V = 500 m/sec) for pumice. 





Figure 30.- Altitude at maximum luminous intensity for pumice. 
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Figure 33 *" Front face peak: temperature, Ti, for iron. 







= 70 km/sec 



Figure 35* - Front face melting zone for iron. 




=70 km/sec 



Figure 37.- Final altitude (V = 500 m/sec) for iron. 



iron 




= 70 km/sec 



Figure 39* - Altitude at maximum luminous intensity for iron. 









